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liquid-gas transition, the dynamical density fluctuations around the CP is analyzed using relativis- 
tic fluid dynamics for a viscous fluid. It is shown that relativistic effects on the dynamical structure 
factor of the density fluctuation appears only in the width of the sound and thermal modes through 
a modification of the transport/kinetic coefficients. It is found that the mechanical density fluctu- 
ation which is enhanced by the relativistic effects are attenuated, whereas the entropy fluctuation 
in turn becomes more prominent around the QCD CP. This attenuation of the density fluctuation 
around the QCD CP may imply that if the suppression or even total disappearance of Mach cone 
is observed as the incident energy of the relativistic heavy-ion collisions is lowered, it can be a 
signal of the existence of the QCD critical point. 



5th International Workshop on Critical Point and Onset of Deconfinement - CPOD 2009, 
June 08 - 12 2009 

Brookhaven National Laboratory, Long Island, New York, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Identifying the QCD Critical Point 



Teiji Kunihiro 



1. Introduction 



An interesting feature of the QCD phase diagram is the existence of a critical point[[T], 0]; 
see Refs. [|], |], ||] for possible variants and alternatives. At the QCD critical point (CP), the first 
order phase transition terminates and turns to a second order phase transition. It is now known 
that the QCD CP belongs to the universality class H in the classification scheme by Hohengerg- 
Halperin[^], i.e., the same as that of the liquid-gas phase transition point[Q, [Sj] _ This implies that 
the density fluctuating mode and, generically, fluid dynamic modes coupled to conserved quantities 
are the softening modes at the CP. 

The a mode that is coupled to the fluctuation of the chiral order parameter (qq)) becomes the 
soft mode of the chiral transition at finite temperature T but vanishing chemical potential pL = 0. 
The QCD CP exists, however, at finite chemical potential /I and for finite current quark masses. In 
such a case, as is first shown in [|9j], the scalar- vector cross term (: (q'f > q){qq) '■) does not vanish, 
and, hence, the scalar mode is coupled with the density-density correlator (: (q^q) 2 :); this is 
because charge conjugation symmetry is violated with finite jj,, and the left- and right-handed 
quarks are coupled owing to the breaking of chiral symmetry. 1 Then, the would-be soft mode, the 
a mode, remains massive (a fast mode) and becomes a slaving mode of the genuine soft mode 
given by the fluid dynamical modes coupled to the density fluctuation^], 0]. Here we note that the 



diverging behavior of the density fluctuation around the QCD CP was first shown in [10]; see also 



[11] 



Recently, it has been argued [12] on the basis of the analysis of the liquid-gas phase transition 
by Onuki[]T3|] that the bulk viscosity may show a singular behavior around the QCD CP. 

We explore how the singularities of the thermodynamic quantities as well as the transport 
coefficients affect the dynamical density fluctuations around the QCD CP thereby we can have a 
hint how to identify the signal of the existence of QCD CP by experiment[[l4|]. Our investigation 
is based on an explicit use of the relativistic fluid dynamic equations for a viscous fluid. We also 
demonstrate that even the so called first-order relativistic fluid dynamic equations have generically 
no problem to describe fluid dynamical phenomena with long wave lengths contrary to a possible 
naive suspect. 

We find that the mechanical density fluctuation which is enhanced by the relativistic effects 
are attenuated, but the entropy fluctuation in turn becomes more prominent around the QCD CP, 
which is, to our surprise, precisely the same in effect as the critical behavior in the nonrelativistic 



case. [ |15| ] On the basis of these findings, we speculate that if a Mach cone, which originates from 
the dynamical density fluctuations stimulated by a fast particle and seems to have been observed 
in RHIC experiment[|T^], disappears or is strongly suppressed as the incident energy of relativistic 
heavy-ion collisions is lowered, it can be a signal of the existence of the QCD CP provided that the 
incident energy is still high enough to make jets[|l7|]. 



'it is also notable that the density-density correlator is a part of the polarization tensor {: {qy^ q){qj v q) :), which 
suggests that the vector coupling g v (qy^q) 2 plays an important role for the static and dynamic properties of the QCD 

cp. [Si 
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2. Dissipative Relativistic Fluid Dynamics 

Relativistic fluid dynamic equations consist of the balance equations for energy-momentum 
and particle number, d^T^ v = and d^N^ = where T^ v is the energy-momentum tensor and N 11 
the particle current: 

r' iV = (e+iV i « v -jy tv + T' iV ) N"=nu^ + V^. (2.1) 

Here e is the energy density, P the pressure, the flow velocity, and n the particle density; the 
dissipative part of the energy-momentum tensor and the particle current are denoted by T^ v and 
v^, respectively. 



The so called first-order equations such as Landau[18] and Eckart[19] ones are parabolic and 
formally violates the causality, and, hence, are called acausal. Moreover, the Eckart equation which 
is defined for the particle frame where the particle current does not have a dissipative part (v^ = 0) 



shows a pathological property that the fluctuations around the thermal equilibrium is uns table [|20|]. 



The causality problem is circumvented in the Israel-Stewart equation||21|], which is a second-order 
equation with relaxation times incorporated. 

One should, however, notice that formally acausal fluid dynamic equations should be valid in 
describing fluid dynamical phenomena with longer wavelengths than the mean free path. In fact, 
we will see that the results for fluid dynamical modes with long wave lengths are qualitatively the 
same irrespective whether the second-order or first-order equations are used. 

As for the instability seen in the Eckart equation, a new first-order equation in the particle 
frame constructed by Tsumura, Kunihiro and Ohnishi (TKO) [|22|] has no such a pathological 



behavior n23|]. We employ Landau[[T^], Eckart[[T9|] and Israel-Stewart(I-S)pT|] equation as typical 
equations, and TKO equation in particle frame. 



3. Dynamic structure factor of density fluctuations 

In this section, we first present a procedure for deriving the dynamical structure factor (spec- 
tral function) of the density fluctuation for the Landau equation. Then we just give the results 
for Eckart, Tsumura-Kunihiro-Ohnishi(TKO) equation and the Israel-Stewart equation in particle 
frame. 



3.1 Landau equation (energy frame) 

The dissipative terms in the Landau equation are given by 



nT 
w 



(3.1) 



where r\ is the shear viscosity, £ the bulk viscosity, K the thermal conductivity and w = e+P 
the enthalpy density. A^ v = g^ v — u^u v is the projection operator on the space-like vector, and 
= A^ v d v the space-like derivative (gradient operator). 

We calculate the dynamical structure factor, or the spectral function, of the density fluctuation 
around the thermal equilibrium state, as was done for the non-relativistic case[24, 25]; we refer to 
the comprehensive text book Ref. [25]. 
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Let us express the dynamical quantities as follows, n(x) = iiq + 8n(x), e(x) = £o + 8e(x), 
P{x) = Po + 8P(x), ;U(x) = Ho + 8jl{x), and u^(x) = + 8u^(x), where the variables with the 
subscript o denote those in the equilibrium state in at rest; u.q = (1,0). Then, 8u tl {x) = (0, 5v(jc)) 
with 8\(x) yet to be determined together with other quantities like Sn etc. Now we have five equa- 
tions for seven unknown quantities, 8n, 8T, 8P, 8s, and 8\. Choosing Sn and ST as independent 
variables, we have SP(x) = ^Sn{x) + ^^8T{x) and Ss(x) = -^§^8n{x) + %8T{x). Here 
Cn = To(ds/dT)„ and dp = To(ds/dT)p are the specific heats at constant density and pressure, re- 
spectively, c s = (dP/de) s the sound velocity, ctp = — {\ /n ){dn/ dT) P the thermal expansivity 
at constant pressure, and y = c P / c„ the ratio of the specific heats. The equation describing the 
density fluctuations is finally found to be 

{ |_ _ K M V 2}5n + n V • 8y+ (1 - C -^)V 2 ST = 0, (3.2) 

W0 ^l _ ^ v 2 5v - (C + irj)V(V • 5v) + ^ V5n + VSr = 0, (3.3) 

a? 3 «o7 7 

vvocjap a K SL v 2 ]Sn + { ^|_ + K{ ±? _ ±) V 2 }Sr = 0. (3.4) 
1 n 7 dt my s Tq dt y Tq 

Now we are interested in the dynamical structure factor of the density fluctuation as given by 

S nn (k,co) = (Sn(k,co)Sn(k,t = 0)), (3.5) 

where 5«(k, w) is the Fourier transform of the density fluctuation, and ( ) denotes the thermal 
average in the equilibrium. 

Note that ST and Sn are statistically independent in fluid systems which is established in Ein- 
stein fluctuation theory [j25l]; (8T(k,0)8n(k,0)) =0, and similarly, (5vn (k,0)5n(k,0)) =0, where 
5v||(k,0) denotes the longitudinal component of the velocity field. Thus we have the dynamical 
structure factor 

Snn(K(o) 1, 2r R £ 2 i r B k 2 r B k 2 



((8n(k,t = o)) 2 } K r co 2 +rlk 4 r {co-c s k) 2 +r 2 k 4 (co+c s k) 2 +r 2 k 4 ~ 

Here, T R = K/(n c P ) and T B = £fc(y- 1) + V/] +jC 2 T ( y ± - 2^a P ) with v, = (£ +4tj/3)/w 
being the (relativistic) longitudinal kinetic viscosity; notice wo in the denominator. We also remark 
that the width Tr is identified with the thermal diffusivity %. 

We see that the spectral function has three peaks at frequencies CO = and CO = ±c s k: The peak 
at CO = corresponds to thermally induced density fluctuations, which is called Rayleigh peak; 
while the two-side peaks at CO = ±c s k correspond to mechanically induced density fluctuation, i.e. 
sound waves. These two peaks are called Brillioun peaks. 

We find that relativistic effects manifest themselves only in the width of the Brillouin peaks 



Tb, whereas the width of the Rayleigh peak is the same as the non-relativistic case[|24j, |25|]. The 
relativistic effects in Tb appear in two ways, as seen by the expression, 

r B = i1f + sr!f (3.7) 

where R = \\x{j— 1) + V/] and 8T\f = \c 2 s Tq(k/wq - 2%ap). Firstly, the first term has a non- 
relativistic counter part only with a replacement of the enthalpy density wq with the mass density po 



4 



Identifying the QCD Critical Point 



Teiji Kunihiro 



300 
250 
200 
150 
100 
50 





Landa 






Minim 


_ 









•0.1 -0.08-0.06-0.04-0.02 0.02 0.04 0.06 0.08 0.1 
to[1/fm] 





1.1 -0.08-0.06-0.04-0.02 0.02 0.04 0.06 0.08 0.1 
to[1/fm] 



Figure 1: Left panel: The solid line is the spectral function for the density fluctuations from Landau/Israel- 
Stewart equation with k = 0.1[l/fm], /i Q = 200 [Me V], T = 200[MeV], f]/(n Q s ) = C/(«o^o) = 0.3 and 
kTq I («o*o) =0.6; the dashed line in the minimal relativistic case. Right panel: the result from TKO equation 
with the same parameters. Notice that the vertical scale is different from that of the left panel. 



in V/. We call this modification the minimal relativistic (MR) effect. On the other hand, the second 
term 8r]f is a genuine relativistic effect that is absent in the non-relativistic case; this comes from 
the mass-energy equivalence inherent in the relativistic theory. 

Figure [j] shows the spectral function Eq.(3i)) and the minimal relativistic case with 8Tb = 



for the parameter set k = 0.1[l/fm], jUo = 200[MeV], T = 200[MeV], T7 /(n s ) = £/{n s ) = 0.3 
and kTo/(«oSo) = 0.6. Note that «o^o represents the entropy density in the equilibrium state because 
so is the entropy per particle number. As is expected, Figjl] shows that the Brillouin peaks owing 
to the sound mode is enhanced by the relativistic effects, while the Rayleigh peak owing to the 
thermal mode is the same as in the non-relativistic case. 

3.2 Stable first-order fluid dynamic equation in the particle frame 

In the case of TKO equation^] in the particle frame (v^ = 0), the dissipative terms are given 

by 

t" v = 7] [d£u v + - ^ V (d ± -u) ] - C'(3^ M V - A^ v )(d ± -u) + K{u^dlT + u v d£), (3.8) 
where £ = £/ (3y — 4) 2 . From the same procedure as taken for the Landau equation, we obtain the 



dynamical structure factor with the same form as (3.6) but with a different form of the width of the 



Brillouin peaks [14]; 



r B = \[%{y- i) + vI KO -^4(KT + 3a = r™. (3.9) 

2 n c P 

We emphasize that since the fluid dynamic fluctuations around the equilibrium state is relaxing 
in TKO equation, we have obtained the dynamical structure factor without any obstruction, in 



contrast to the case of Eckart equation[19], for which we have found[14] that the unstable behavior 
of the density fluctuation around the thermal equilibrium prevents us from obtaining the dynamical 
structure factor in a sensible way. 
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We see that the genuine relativistic effect in the width of the Brillouin peaks reads dT^ = 
— ( kTq + 3 £' ) , which is definitely negative, implying that the relativistic effect acts to enhance 
and sharpen the spectral function of the density fluctuation in comparison with the minimal rela- 
tivistic case than in the energy frame; see the right panel in Figjl]; notice that the scales of the right 
and left panel are different. 

3.3 Israel-Stewart equation 

The derivation of the dynamical structure factor for the Israel-Stewart (I-S) equation[|T]] can 
be performed much the same way as for the Landau equation with some additional complications 
due to the presence of the relaxation times. However, it is found that if the relaxation time j8i of 
the heat current is so small as to satisfy the condition /3i < the I-S equation in particle frame 
takes over the pathological behavior of Eckart one even though the relaxation time is finite; i.e., the 
density fluctuation will not relax down. 

If we assume that the relaxation time is sufficiently large to satisfy the inequality j3i > then 
we can have the spectral function, 

1 ?vk 2 1 Vvk 2 
S nn (k,co)/((8n(k,t = 0)) 2 )= (1--) , * +-[ 



Y'(0 2 + x 2 k 4 J ((0-c s k) 2 + r 2 k 4 



+ 0{k z ) x [- 



(co+c s k) 2 +r 2 k^ " ' L fl)2 + i/(A)C) 2 

l/fcq 2wo/[(j3iwo-l)K:ro] 
^co 2 +\/(2^) 2 co 2 + W y[{p iWQ -l)KT f- K - ) 

Here, j8o and fi 2 denote the relaxation times for the bulk viscosity and the shear viscosity, respec- 
tively. Apparently the spectral function has six peaks including the conventional three peaks, but 
the new three Lorentzian functions should vanish in the long wavelength limit k — > 0, because the 
strength of these is in the second order of k. Therefore Israel-Stewart equation gives completely the 
same result for the dynamical structure factor for the density fluctuations in the long wavelength 
limit as that Landau equation gives. That is, the relaxation times do not affect the result in the fluid 
dynamical regime. 



4. The behavior around the QCD critical point 

We are now in a position to analyze the behavior of the spectral function of the density fluctu- 
ations around the QCD critical point, on the basis of the dynamic as well as static scaling laws for 
the liquid-gas transition [15, [L3| ] . 



The specific heat at constant density c n and the the isothermal compressibility Kj = ^-{jp) T 
diverge at the critical point as c„ = cot~ a and Kj = Kq(~^, where t = \(T — T c )/T c \ is a reduced 
temperature. The values of these critical exponents are known to be a ~ 0.11 and y~ 1.2, respec- 
tively. 

At the critical point, the correlation length £ diverges as E, = ~ v with v ~ 0.63. The bulk 



viscosity may also show a singular behavior around the QCD CP, as shown in Q13| , |12| ] ; £ = at > 
with £o being a constant. The exponent for the liquid-gas transition is predicted [ |l3j ] to be zv — a 
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with the dynamical critical exponent z; z ~ 3. Thus we have a? ~ 1.8. The singular behavior of the 
thermal conductivity around the critical point is given by a K ; K = Kot~" x with a K ~ 0.63 where Kq 
is a constant. 



Now using these formulae, we can showQ14[] that the width of the Rayleigh peak behaves as 
r R ~ t^" K , which tells us that the width Tr becomes narrow as the QCD CP is approached. We 
emphasize that this result is independent of the choice of the relativistic fluid dynamic equation or 
frame. 

We find that the width of the Brillouin peaks Tb show the following critical behavior for the 
Landau and TS equations, Tb ~ J^t ~ a <>. We note that this singularity comes from that of the bulk 
viscosity. In the case of TKO equation, we first note that the critical behavior of the effective 
bulk viscosity is given as t, = p^zjp ~ t 2 Y~ a z, which shows that the effective bulk viscosity 
has a positive exponent and does not show a singular behavior because ar ~ 1.8 and y ~ 1.2. 
Instead, the singularity of the Brillouin peaks for TKO equation comes from that of the thermal 
conductivity [Q]; T B ~ H a *-«). 

Anyway, we have found that the width Tb diverges at the QCD CP, irrespective of the rela- 
tivistic fluid dynamic equations, although the strength of the singularity may differ depending on 
the choice of the fluid dynamic equationJH]]. 
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Figure 2: Left panel: The spectral function at t = \T — T c \/Tc = 0.4 and k — 0.1 [1/fm]. The solid line 
represents the Landau and I-S cases. The dashed line represents the TKO case. The strength of the Brillouin 
peaks becomes small due to the singularity of the ratio of specific heats. Right panehThe spectral function 
at t = 0.1 and k — 0.1 [1/fm]. We see that the Brillouin peaks which correspond to sound wave dies out and 
the difference between the Landau and TKO cases disappears. 



It turns out that these singular behaviors of the width of the Brillouin peaks around the QCD 
CP may not be observed: Note that the strengths of the Rayleigh and the Brillouin peaks are given 
in terms of y, which behaves like y = c p /c n ~ t~7 +a — > oo, in the critical region. Then the strength 
of the Brillouin peaks is attenuated and only the Rayleigh peak stands out in the critical region, as 
follows; 



S„„(k,G))~((Sra(k, 0)) 2 > 



2T R k 2 
(0 2 + r 2 R k 4 ' 



(T 



(4.1) 



Figure |2] shows how the dynamical structure factor for the density fluctuations in the Landau (I- 
S) and TKO cases behaves around the QCD critical point for k = 0.1 [1/fm], t = 0.5 and t = 0.1, 
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respectively. We see that the strength of the Brillouin peaks becomes small and dies out as the 
system approach the QCD CP. In addition, the static correlation function ((8n(k,t = 0)) 2 ) shows a 
singular behavior in the forward direction k = 0. This is known as the critical opalescence [|l5|, |25| ]. 
Then the strength of the Rayleigh peak will be most drastically enhanced in the forward angle in 
the critical region. 

Why at all do sound modes or mechanical density fluctuations die out at the critical point? To 
answer this question, we first note that the fluid dynamic regime is expressed as t, « X s , where 
X s represents the wave length of the sound mode. When this condition is satisfied, the sound mode 
can be developed^]. However, in the vicinity of the critical point, the correlation length E, becomes 
very large and eventually goes to infinity, so the above inequality can not be satisfied, and the sound 
mode is hardly developed in the vicinity of the critical point[^]. This is the reason why the sound 
mode or density fluctuation is hardly developed around the QCD critical point. 



5. Possible disappearance or suppression of Mach cone around the QCD CP 

The result in the last section suggests that phenomena inherently related to the existence of the 
mechanical density fluctuations may disappear around the critical point. One of such phenomena 



is the possible Mach cone formation [17] by the particle passing through the medium with a speed 
larger than the sound velocity c s . Such a Mach-cone like particle correlations are observed in the 
RHIC experiment[jl^]. If such three-particle correlations have been confirmed to be a Mach-cone 
formation, then the disappearance or suppression of the Mach cone would be a signal that the 
created matter has passed through the critical region, showing the existence of the QCD critical 
point. Even if the thermal wake also contributes to the formation of Mach cone, a suppression of 
Mach cone may be expected by the attenuation of the dynamical density fluctuations. So it would 
be very interesting to see possible variation of the strength of the Mach cone according to the 
variation of the incident energy of the heavy-ion collisions. In theoretical side, it is an intriguing 
task to explore the fate of Mach cone with an equation of state that admits the existence of the CR 



6. Summary and concluding remarks 

Motivated by the fact that QCD critical point (CP) belongs to the same universality class as 
the liquid-gas transition, we have analyzed the dynamical density fluctuations using relativistic 
fluid dynamics. We have found that the mechanical density fluctuations are attenuated owing to 
the divergence of the correlation length ^ around the QCD CP; on the other hand, the entropy 
fluctuation in turn gets enhanced and tends to makes a single peak around the QCD CP in the 
dynamical structure factor of the density fluctuations. Sunch an attenuation of the mechanical 
density fluctuations may lead to a suppression or even total disappearance of Mach cone at the 
QCD CP. If the Mach cone formation is confirmed at the incident energy ^/sm = 200 GeV in 
RHIC experiments, possible disappearance or strong suppression of a Mach cone along with the 
lowering of the incident energy can be a signal of the existence of the critical point; the created 
matter should have gone through the critical region of the CP. 
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Eexplicit calculations with equation of motion which admittes the existence of the critical 
point is necessary for confirm the fate of Mach cone formation. To make a direct connection with 
RHIC experiment, we should analyze the density fluctuations with the expanding back ground[26]. 

As a future work, we should study the coupling between the thermal fluctuations and the 
transverse mode using the mode-mode coupling theory in the close vicinity of the critical point for 
the relativistic case. 

Finally, we would like to indicate that there are still other possibilities in the structure of the 
QCD phase diagram [|, ||] : For example, there may exist multiple critical points when the color 
superconductivity and the vector interaction like g v (qY^q) 2 [[J and/or the anomaly term[Q] are 
incorporated. It suggests that the QCD matter is very soft along the critical line when the color 
superconductivity is incorporated, which may imply that there are large fluctuations of various 
observables including diquark-density mixed ones. 

T.K. is grateful to Krishna Rajagopal for his valuable comments. This work was partially 
supported by a Grant-in-Aid for Scientific Research by the Ministry of Education, Culture, Sports, 
Science and Technology (MEXT) of Japan (Nos. 20540265, 19-07797), by Yukawa International 
Program for Quark-Hadron Sciences, and by the Grant-in-Aid for the global COE program " The 
Next Generation of Physics, Spun from Universality and Emergence " from MEXT. 
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